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ABSTRACT 

We observed the globular cluster NGC 6652 with Chandra for 47.5 ks, detecting six known X-ray 
sources, as well as five previously undetected X-ray sources. Source A (XB 1832-330) is a well-known 
bright low-mass X-ray binary (LXMB). The second brightest source, B, has a spectrum that fits 
well to either a power-law model (F ~ 1.3) or an absorbed hot gas emission model (kT ~ 34 keV). 
Its unabsorbed 0.5-10 keV luminosity {Lx = 1.6~^ 1 x 10 34 erg s _1 ) is suggestive of a neutron star 
primary; however, Source B exhibits unusual variability for a LMXB, varying by over an order of 
magnitude on timescales of ~ 100 s. Source C's spectrum contains a strong low-energy component 
below ~ 1 keV. Its spectrum is well fit to a simplified magnetic cataclysmic variable (CV) model, 
thus the soft component may be explained by a hot polar cap of a magnetic CV. Source D has an 
average Lx (0.5-10 keV) ~ 9xl0 32 erg s _1 and its spectrum is well fit to a neutron star atmosphere 
model. This is indicative of a quiescent neutron star LXMB, suggesting Source D may be the third 
known LMXB in NGC 6652. Source E has L x (0.5-10 keV) - 3xl0 32 erg s"\ while Source F has 
Lx (0.5-10 keV) ~ lxlO 32 erg s _1 . Their relatively hard X-ray spectra are well-fit by power-law or 
plasma emission models. Five newly detected fainter sources have luminosities between 1-5 x 10 31 erg 
s . NGC 6652 has an unusually flat X-ray luminosity function compared to other globular clusters, 
which may be connected to its extremely high central density. 

Subject headings: binaries : X-rays — cataclysmic variables — globular clusters: individual (NGC 
6652) — stars: neutron 



1. INTRODUCTION 

Low-luminosity X-ray sour ces in globular c l usters 
were identified with Einstein (jHertz fc Grindlavl fl983l ) 
as a separate class of X-ray sources from bright 
(10 36 <Lx<10 38 ergs/s) LMXBs containing neutron 
stars (NSs) . They suggested that low-luminosity sources 
were primarily cataclysmic variables (CVs), with simi- 
lar accretion rates to the LMXBs but potential energy 
wells 1000 times shallower, but al so including some NS 
LMXBs in quiescence (qLMXBs). IVerbunt et all (fl98l 
argued that the brightest low-Lx systems (Lx> 10 33 
ergs/s) were too bright to be CVs, and must be qLMXBs. 
The answer, from Chandra observations, has been mixed; 
some of the brighter \ow-Lx systems have been identi- 
fied with qLMXBs, an d some with CVs (e.g. in 47 Tuc, 
IGrindlav et "aTJ l205lah . 

The brighter low-Lx systems (5 x 10 32 <L.y<10 35 
ergs/s) are of special interest. The bright CVs in this 
range are at the top of the X-ray luminosity distri- 
bution of CVs. As disk-fed CVs transition to opti- 
cally thick fl ows at high mass transfer rate s, suppressing 
hard X-rays ([Patterson fc Raymondlll985h . it is thought 
that reaching high X-ray luminosities may require mag- 
netically channeled accretion (permittin g optically thin 
flows) onto massive whi te dwarfs fe.g. IGrindlav et al.l 
119951 : llvanova et al.|[2006() . The relative lack of dwarf no- 
vae outbursts in globular clusters als o supports the map 
netically channeled accretion idea (jShara et all I199€ 
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However, there is little direct evidence for magnetic CVs 
in globular clusters; published evidence consists of A4686 
He II emission (associated with intermediate polar CVs, 
in which the accretion disk is trun cated by the primary 's 
B field) from 3 CVs in NGC 6397 (Edmonds et al.l[l999h . 
and very so ft blackbody-like co mponents in the CVs X10 
in 4 7 Tuc dHeinke et all 120(351) and IE 1339.8+2837 in 
M3 (|Stacev et al.ll2011[ ). Further evidence of the nature 
of the most luminous CVs, in particular direct evidence 
of strong magnetic fields, will be important. 

The most luminous qLMXBs (10 35 > L x > 10 34 
ergs/s) are brighter than t he expected emission from 
residual heat in the crust (jYakovlev fe Pethickl 120041 ) , 
so must be powered by continuing accretion. However, 
they are fainter than expected for X-ray binaries under- 
going outb ursts driven by the standard disk instability 
model (e.g. 11^120001: IWiinands et al.ll2006l) . We don't 
yet understand the behavior of LMXBs that maintain 
these intermediate luminosities for years, s uch as XMMU 
J174716. 1-281048 ([DeTSanto et al.l 1200% and M15 X-3 
(|Heinke et al.l 120091 ). The neutron star's magnetic field 
is expected to exercise a "propeller effect" which should 
stop the accretion of material onto the neutron star sur- 
face at low mass transfer rates (|Illarionov fc Sunvaevl 
Il975f ). On the other hand, simulations of the interaction 
of weak accretion with a magnetic prop eller predict some 
mate rial will reach the neutron star (e.g. lRomanova et al.l 
120021 ). suggesting that NS LMXBs in a propeller state 
might produce accretion at intermediate luminosities. 

NGC 6652 contains a moderately bright (L x = 10 36 
ergs/s) LMXB (A). Three other sources were identi- 
fied in a 1.6 ksec Chandra HRC observation in 2001 
(|Heinke et all 1200 If) . The brightest faint source, B, at 
(0.5-7)~ 6 x 10 33 ergs/s, has an optically (Hubble 
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Spac e Telescope, HST) id entified blue variable counter- 
part (jDeutsch et al.lll~998l ). It is UV-bright but appears 
on the main sequence in a V vs. V — I color-magnitude 
diagram, suggesting a quiescent LMXB with a weak disk 
(jHeinke et all 120011 ). A 5-ks 2008 Chandra observation 
showed rapid variability, from 2 x 10 33 up to > 5 x 10 34 
crgs/s on timescales <2 minutes, not previously seen 
from LMXBs at such low luminosities. Spectral fitting 
of countrate-selected bins suggested spectral variability, 
but it was unclear whether the behavior was true flaring 
or variations in obscur ing column from an edge-on source 
(ICoomber et al.l[20Tl . 

The next brightest source, C, is very blue in V — 
I HST observations, which combined with its Lx of 
6 x 10 32 ergs/s, indicat es it is likely a very luminous 
CV (jHeinke et all 120011) . The 5-ks Chandra observa- 
tion found an extremely so ft spectrum, with mo st of the 
90 counts below 0.8 keV (jCoomber et all 120111) . D, at 
Lx = 5 x 10 32 ergs/s in the (optically crowded) cluster 
core, also showed a soft spectrum, with a suggestion of a 
line in the (73-count) spectrum. The 5-ks Chandra ob- 
servation revealed three other faint sources, of which G 
lies outside our field of view. 

NGC 6652 has been con sidered to be a cluster of only 
moderate ce ntral density (lHarrisI [20To() and stellar en- 
counter rate (|Verbuntll2003t ). making its possession of a 
bright LMXB somewhat unu sual. However, a recent sur - 
face brightness profile from iNovola fc Gebhardtl ()2006l ) 
finds a sharp core with a significantly higher central sur- 
face brightness than previously measured, identifying it 
as a core-collapsed cluster. We have computed the cen- 
tral luminosity densities of Galactic globular clusters us- 
ing the Noyola & Gcbhardt parameter updates, and note 
that NGC 6652's central density is the fifth highest. 

2. DATA REDUCTION 

We observed the globular cluster NGC 6652 on June 
3rd, 2011 with the Chandra ACIS-S detector, using two 
active CCDs (S2 and S3) with a 1/8 subarray to abate 
pilcup from Source B. The total observation time was 
47.5 ks. The data were reduced using CIAO version 
4.S0. A new bad pixel file was generated with the CIAO 
acis ..runJiotpix tool. The observations were further pro- 
cessed using acis -process -events and filterecO in order to 
create level 2 event files, as detailed in the CIAO Science 
Threads^. 

We ran the CIAO wavdetect toolQ over an energy range 
of 0.5-10 keV. From this, we found six X-ray sources, all 
located within the cluster half-mas s radius (Figure d!). A 
seventh previously detected source (|Coomber et al.l2~011[ ) 
lay outside our field of view. The source positions agree 
with our previous 2008 Chandra observation. In addi- 
tion, we ran wavdetect over several other ranges (0.3-7 
keV, 0.3-2 keV, 2-7 keV); however, no previously undc- 
tected sources were fo und. We further ran PWDetect 
(|Damiani et al. 1997a. b) over a range of 0.3-7 keV. Both 
PWDetect and wavdetect use wavelet transforms to lo- 
cate possible sources; we have found that wavdetect ef- 

4 http:/ /cxc. cfa.harvard.edu/ciao/ 

5 http:/ /cxc. cfa.harvard.edu/ciao/thrcads/creatcL2/ 

6 http:/ /cxc. cfa.harvard.edu/ciao/threads/all. html 

7 http:/ /cxc. harvard.edu/ciao/thrcads/wavdetcct/ 

8 http://www.astropa.unipa.it/progetti_ricerca/PWDetcct/ 



■ 

■ / 


\ 






„(§) 






















48.0 


47.0 46.0 18:35:45.0 44.0 
Flight ascension 


43.0 



Fig. 1— ACIS-S image of NGC 6652 over 0.3-7 keV, with X-ray 
sources are labeled. The core radi us (0.10' ) and half-mass radius 
(0.48' ) are shown ((Harris 2010). A previously detected X-ray 
source (Source G) lies outside our field of view (see ICoomber et al.l 
1120111 )). while Source L lies ~1.4' from the center of the cluster. 

ficiently finds sources without spurious detections over 
large uncrowded fields of view, while PWDetect can 
more effectively identify faint sources near bright sources, 
though it has a higher spurious de tection rate, making 
it less suitable for lar ge fields (e.g. iHeinke et al.l l2003"at 
iBogdanov et al.l l2010T l . In order to reduce the num- 
ber of spurious detections, we only include detections 
above 5 a. We identify four previously undetected X- 
ray sources within the cluster half-mass radius (Sources 
H through K) and one source outside the cluster half- 
mass radius (Source L). Our detection limit of 10 counts 
translates to roughly 1 x 10 31 ergs/s (0.5-2.5 keV), or 
2 x 10 31 ergs/s (0.5-10 keV) for standard absorbed pow- 
erlaw spe ctra with ph oton indcx=2 at the 10 kpc cluster 
distance (|Harrisll20l6T ). 

We selected events within a circle around the sources 
(radius 1 " for H, J, K and L, radius 0.75 " for I, due to its 
proximity to the readout streak), and background from 
local regions, determining the background subtracted 
source counts. The positions and net counts for each 
observed source are listed in Table [TJ with 95% confi- 
dence posit i onal u ncertainty values from equation (5) of 
iHong et al.l (|2005[ ). For Sources B through F, we selected 
events from a circular source region (2" for B, 1.4" for 
C through F) centered on the wavdetect positions. For 
Sources B, E and F, an annulus of inner radius 5" that 
excluded the source was chosen for the background. Due 
to the proximity of Sources C and D, a larger annulus 
of inner radius 6 " that excluded both sources was cho- 
sen for the background. Using these regions, the CIAO 
specextract script was run to extract source and back- 
ground spectra, generating RMF and ARF files. XSPEC 
version 12.70 was used to examine and fit the spectra 

9 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspcc/ 
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Fig. 2. — ACIS lightcurve extracted from Source B in the energy 
band 0.3—7 keV, using a binning of 100 s. The lightcurve is split 
into three consecutive periods (top, middle and bottom) in order 
to clearly display the peaks in the lightcurve. A corresponding 
variability index of 10 (variable) was found using the CIAO glvary 
tool. 

(discussed in the following sections). 

3. X-RAY ANALYSIS 

3.1. Source A 

We extracted 1537 counts from Source A's readout 
streak (using a 99.5 by 3.5 pixel box region), group- 
ing them into 14 bins with a minimum of 100 counts 
per bin. We were able to successfully (x 2 =l-05 for 11 
degrees of freedom) fit Source A's spectrum to an ab- 
sorbed power-law model (r = 1.9 ± 0.3), with Lx(0.5-10 
keV)=7.0to;7 x 10 35 er g s / s - A hydrogen column density 
of 27i" x 10 20 cm 2 was measured, significantly larger 
than the cluster value of 5xl0 20 cm" 2 (jHarrisI [20Toh . 
Previous studies of A's X-ray spectrum with ASCA, Bep- 
poSA X. and XMM (iM^kai 'fc Smald l2000HParmar et al.l 
[20011: ISidoli et al.ll2008D agree in requiring additional in- 
trinsic absorption, with their single-component best-fit 
spectra resembling ours. Those previous observations 
(optimized for this bright source) were able to better 
constrain A's spectrum, measuring partial covering and 
multiple spectral components, which our data are insuf- 
ficient to constrain. 

3.2. Source B 
3.2.1. Timing Analysis 

A 474 bin lightcurve (100 s/bin) was extracted for 
Source B using the ClAO dmextract tool (Fig. [5]). 
The observation shows clear variability by factors over 
an order of magnitude on timcscales of ~ 100 s, un- 
characteristic of low-mass X-ray binaries. Wc verified 
the variability with the CIAO glvary tool, which ap- 
plies the Grcgory-Loredo algorithm. The algorithm sep- 
arates the observation into bins based upon time and 



look s for significant variation , returning a variability in- 
dex ()Gregorv fc Loredolll992t) . For Source B, a variabil- 
ity index of 10 was found, corresponding to a variabil- 
ity probability of ~ 1.0. Power spectra (generated us- 
ing XR0N0S3) show no clear periodicity in Source B's 
X-ray emission (using either 100 s or 0.54 s binning). 
In order to test whether the variation may be due to 
changes in the obscuring column, the observation was 
filtered into eight separate ranges based upon the count 
rate. The count rate ranges were chosen to ensure each 
range contains at least 400 counts. The hardness ratio 
within each range was then quantified by examining the 
ratio of the number of counts with photon energy be- 
tween 0.3-1 keV to the number of counts with photon 
energy between 0.3-7 keV (see Table [2]) . If the varia- 
tion was due to obscuration, we would expect to find 
the lowest hardness ratios in the lowest count rate range. 
This is seen to some degree as the lowest hardness ratio 
(0.16±0.02) is found between 0.05-0.1 counts s _1 and the 
largest hardness ratio (0.23±0.02) is seen between 0.27- 
0.35 counts s . Wc further examined the spectrum for 
each of these ranges (below). 

3.2.2. X-ray Spectral Analysis 

Source B's total spectrum was binned (using the GRP- 
PHA tool) so that each bin contained a minimum of 100 
counts. We included a photoelectric absorption compo- 
nent (XSPEC model phabs) with each model, leaving the 
hydrogen column denstiy, Nu, as a free parameter. The 
metal abundance for the mekal model was kept at [Fc/H] 
= -0.81 jHarrH EIHSEE Considering a range from 0.5 
to 10 keV, Source B's spectrum was fit reasonably well 
to both an absorbed power-law model (XSPEC model 
pegpwrlw, photon index T = 1.3 ± 0.1) and an absorbed 
hot X-ray plasma model (XSPEC model mekal, tempera- 
ture kT=34l 2 g keV) . We found an unabsorbed X-ray flux 
(0.5-10 keV) of (1.4±0.1) x 10" 12 erg cm" 2 s" 1 from the 
power-law model and (1.3±0.1) x 10~ 12 erg cm" 2 s _1 
from the plasma emis sion model. Using a distance of 
10.0 kpc (|Harrislll996|) . these correspond to luminosities 
of (1.7±0.1) x 10 34 erg s" 1 and (1.6±0.1) x 10 34 erg s" 1 , 
consistent with the 2008 Chandra ACIS-S observation. 
The full results are listed in Table [3J 

We next filtered the spectrum into eight separate 
ranges based upon count rate in order to look for possible 
changes in the spectra. Using a minimum of 30 counts 
per bin, we fit the spectra to an absorbed mekal model 
over 0.5 to 10 keV (FigureEJ. In addition, the five spec- 
tra with the largest count rates contained a pileu p com- 
pone nt (grade morphing parameter fixed at 0.5, iDavisI 
12001 . TheX -ray plasma emission temperature and nor- 
malization, as well as the hydrogen column density (Nh), 
were allowed to vary. The parameters from these fits are 
listed in Table [2j There appears to be a tendency for the 
higher count rate spectra to have lower Nh values, as the 
highest Nh observed appears in the 0.05-0.075 counts s _1 
range and the lowest appears in the 0.27-0.35 counts s _1 
range. However, this tendency is not strong; the maxi- 
mum fitted Nh is only 2.7+J'l x 10 21 cm" 2 , insufficient to 
significantly decrease the countrate. Therefore, we rule 

10 http:/ /heasarc. nasa.gov/xanadu/xronos/xronos. html 

11 Updated 2010; http://www.physics.mcmaster.ca/~harris/mwgc.dat 
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Fig. 3. — X-ray spectra of Source B, with a minimum of 30 
counts per bin. Spectra are filtered by count rate ranges, with best- 
fit absorbed mekal models fit from 0.5 to fO keV. Eight separate 
ranges were fit to produce the results in Table [2] but for clarity 
only four are shown here. From highest to lowest count rate, these 
are >0.35 counts s~ 1 (blue, squares), 0.15 - 0.2 counts s _1 (green, 
stars), 0.075 - 0.1 counts s — 1 (red, triangles) and <0.05 counts s — 1 
(black, circles). The ratio of the data to the model is shown below. 

out that the flux variations are entirely due to obscura- 
tion by changing Nh columns, e.g. due to an edge-on 
accretion disk. 

We placed Source B (as well as the other observed 
sources) on an X-ray color-magnitude diagram (Fig- 
ure |4]). We defined color to be 2.5 log[(0.5-1.5 keV 
counts) /(l. 5-7.0 keV counts)] and plotted it against the 
0.5-10 keV luminosity of the sources. In addition, we 
plotted the X-ray color predictions for power-law and 
mekal models, as well as the X-ray color-luminosity re- 
lation for the nsatmos model. 

3.3. Source C 

A lightcurve (500 s/bin) was extracted for Source C 
(see Figure [SJ- Using the CIAO glvary tool (discussed 
in Section [3.2.ip a variability index of 7 was found, with 
a corresponding variability probability of 0.995. Power 
spectra show no clear evidence of periodicity. The to- 
tal unabsorbed luminosity from the simplified polar CV 
model (below) is L x (0.5-10 keV) = 1.6t°i x 10 33 erg 
s _1 , which is comparable to the 2008 ACIS-S observation 
which found Lx (0-5 -10 keV) = 1.2±%\ x 10 33 erg s" 1 
(jCoomber et al.H201lD . 

Grouping Source C's spectrum into bins with a mini- 
mum of 30 counts each, a notable low energy component 
below ~ 1 keV was found (Figure [6]). At energies below 
~ 0.5 keV, the accuracy of the Chandra response files and 
the sensitivity of the ACIS-S detector are no longer op- 
timal. For this reason, we only considered a range from 
0.5 to 10 keV for our primary fits. Even so, we were 
unable to adequately model Source C to any single com- 
ponent model (see Table[3|). Allowing the Nh to vary, we 
successfully fit the spectrum to a model containing two 
hot X-ray plasma emission components (XSPEC model 
mekal) with temperatures >17 keV and 0.18±0.03 keV. 
Double mekal models are commonl y used to describe the 
spectra of CVs (|Baskill et al.ll2005l) ; however, Source C's 
lower-energy component is unusually strong. In order 
to examine the lower energy emission, we fit the same 
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Fig. 4. — X-ray color-magnitude diagram for detected sources. 
Color is determined as a function of low-energy to high energy 
counts. X-ray luminosity (0.5-10 keV) is plotted a l ong t he ver- 
tical axis. Errors are derived for 1 a from Gchrels (1986). Val- 
ues for Sources A and G are taken from [Coombcr ct al.l l|2011l ) . 
Additionally, the theoretical colors for absorbed power-law, mekal 
and nsatmos models (with the corresponding Lx for nsatmos) are 
shown. 
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Fig. 5. — ACIS lightcurve extracted from Source C over 0.3— 
7 keV, using a binning of 500 s. The CIAO glvary tool gave a 
variability index of 7 (variable). 

double mekal model over 0.2 to 10 keV; however, this 
failed to produce an acceptable fit (xt = 1-99). Source 
C's spectrum was adequately fit to a simplified polar 
CV model, containing a high-tcmpcrature mekal compo- 
nent and a soft blackbody component (e.g. iRamsav et al.l 
I2004D . Again, allowing the Nh to vary, the best fit model 
resulted in a mekal temperature >18 keV and a black- 
body temperature of 70±20 eV. Fitting down to 0.2 keV 
with this model gave a better fit than the double mekal 
model in the same range (jq = 0.95, mekal temperature 
> 14 keV, blackbody temperature of 80±10 cV). 

3.4. Source D 

The lightcurve from Source D showed no evidence of 
variation. A variability index (discussed in Section l5.2.1[) 
of was found, resulting in a low variability probability 
of 0.033. 
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Fig. 6. — Chandra spectrum of Source C fit to an absorbed 
mekal plus blackbody model from 0.5 to 10 keV, with a minimum 
of 30 counts per bin. The dotted lines indicate the two components 
of the model (mekal being of higher energy). The ratio of the data 
to the model is shown below. Resulting parameters from the fit 
are listed in Table [3] 




Fig. 7. — Chandra spectrum of Source D fit to an absorbed neu- 
tron star atmosphere model from 0.5 to 10 keV, with a minimum 
of 30 counts per bin. The ratio of the data to the model is shown 
below. Resulting parameters from the fit are listed in Table [3] 

Source D's spectrum was extracted and grouped to en- 
sure a minimum of 30 counts per bin (Figure [7]). Includ- 
ing a photon absorption component and allowing the Ng 
to vary, Source D was fit to a power-law model resulting 
in a photon index of 5.8t°;2 (see Table©. The fitted N H 
of (42tig) xlO 20 cm" 2 is much larger than the cluster 
value of 5 xlO 20 cm" 2 , and the high power- law photon 
index suggests testing thermal emission models. A mekal 
model proved to be a poor fit; however, a neutron star 
atmosphere model (XSPEC model nsatmos, 1.4 M Q , ra- 
dius of 10 km) proved to sufficiently fit the spectrum. 
From this we found an effective temperature of 111±2 
cV. We used the XSPEC cflux model to find L x (0.5- 
10 keV) = 9.0^:? x 10 32 erg s _1 , consistent with the 
2000 Chandra HRC observation and the 2008 Chandra 
ACIS observation. O ur previous Chandra observation 
(|Coomber et al.|[2~011[ ) suggested the existence of an emis- 
sion line near 1 keV; however, our higher-quality spec- 
trum shows no evidence for this. 



Energy (keV) 

Fig. 8. — Chandra spectrum of Source E fit to an absorbed 
power-law model over 0.5 to 10 keV, with a minimum of 10 counts 
per bin. The ratio of the data to the model is shown below. Re- 
sulting parameters from the fit are listed in Table |3] 

3.5. Sources E + F 

We detected 99 counts from Source E and 35 counts 
from Source F. The relatively low number of counts re- 
duces the accuracy of timing analysis; however, it is still 
possible to use the CIAO glvary tool to look for strong 
variation. A variability probability of 0.78 was found for 
Source E (variability index of 2) and a variability prob- 
ability of 0.45 was found for Source F (variability index 
of 0), thus we do not have strong evidence of variability 
in either. 

The spectrum for Source E was extracted and grouped 
to ensure a minimum of 10 counts per bin (see Figure ©. 
To ensure optimal accuracy from the Chandra response 
files, we considered an energy range from 0.5 to 10 keV. 
Allowing the Nh to vary, Source E was successfully fit to 
both an absorbed power-law model with photon index of 
2-2lg'g and an absorbed mekal model with temperature 
>1.8 keV (Table©. An unabsorbed 0.5-10 keV lumi- 
nosity of 2.9t° 9 7 x 1032 

erg s 1 was determined from the 
mekal model, consistent with our previous 2008 Chandra 
observation. 

The spectrum for Source F was also grouped to ensure 
a minimum of 10 counts per bin, resulting in a total of 
four bins (Figure©. Due to the low number of bins, the 
entire 0-10 keV range was considered. We were unable to 
find an ideal single-component model with which to fit 
Source F. An absorbed nsatmos model proved a poor fit, 
while the null hypothesis probability (NHP) for both an 



absorbed power-law (r = UI^, NHP = 0.31) and an 
absorbed mekal model (temperature > 1.9 keV, NHP = 
0.29) remained high. Using the C-statistic (|Cashlll979l h 
we further examined the suitability of each model by fit- 
ting the unbinned data. For the power-law model 99.9% 
of simulated spectra resulted in lower C-statistic values, 
while 99.7% did for the mekal model and 100% did for the 
nsatmos model (for 10 simulated spectra), suggesting a 
poor fit in each case. The power-law and mekal param- 
eters for the unbinned fitting agree within error to the 



binned fitting parameters (r < 6.0, kT = 1 
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The 



unbinned nsatmos model results in a marginally lower 
temperature of 85ii° eV. The luminosity obtained for 
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Fig. 9. — Chandra spectrum of Source F. The model was deter- 
mined using the C-stastistic, fitting an absorbed mekal model to 
the unbinned data over to f keV. For clarity only, the spectrum 
is shown with a minimum of 5 counts per bin. The ratio of the 
data to the model is shown below. Resulting parameters from the 
fit are listed in Table [3] 

Source F (Table |3|) was consistent with the previous 2008 
Chandra observation. 

3.6. Faint Sources 

The five faintest sources have too few counts for ac- 
curate timing or spectral analysis. Using the number 
of source counts found for each source and the Chandra 
PIMMS version 4.2, we estimated unabsorbed X-ray lu- 
minosities (0.5-10 keV) using a photon index of 1.4 (Ta- 
ble [TJ . There exists marginal evidence for a sixth faint 
source located between C and D; however, it was not 
identified by either wavdetect or PWDetect (possibly due 
to its proximity to C and D). It would be difficult for 
more than one or two X-ray sources over 10 counts (trans- 
lating to Lx(0.5-2.5 keV)> 10 31 ergs/s) to be missed, 
apart from near source A or within 2" of source D. 

3.7. Radial Distribution 

We have analyzed the radial distribution of the cluster 
sources by fitting generalized King mod el profiles using 
maximum likelihood t echni q ues (see, e.g. Grindlav et al.l 
(|200l : iHehTke et all (pOOfih : ICohn et all pOlOh l. In this 
approach, the projected radial distribution is given by 
the expression, 



S(r) = So 



v/2 



(1) 



where a is the large-r power-law slope and ro is a ra- 
dial scale factor that is related to the core radius by r c = 
(2~ 2 / Q — l) 1 ' 2 ro- We assume that the turnoff mass stellar 
distribution, which dominates the optical surface bright- 
ness profile, is described by a King model, i.e. by Eqn. [1] 
with a = —2 and r c = r$. In thermal equilibrium, the 
distribution of the Chandra sources is expected to be de- 
scribed by Eqn.Q]with ax = —3q+l, where q = Mx/M* 
is the ratio of the characteristic source mass to the mass 
of the stars that dominate the optical surface brightness 
profile. We he re adopted the values of r c = 5.8" and 
r h = 29" from iMcLaughlin k van der Marell (|2005f) for 



r (arcsec) 

Fig. 10. — Cumulative radial distribution of Chandra sources in 
NGC 6652 (solid line) with a generalized King model fit (dashed 
line). The sample includes the nine sources within the half-mass 
radius. 

the optical surface brightness profile. It is first necessary 
to correct the source di stribution for b a ckgro und contam- 
ination. Based on the lGiacconi et ail (|2001[ ) extragalac- 
tic source counts, we estimate a total of 2.1 background 
sources within the 1' x 8' subarray. We note that there 
is one source (L) that lies well outside of the half-mass 
radius, which is roughly consistent with this predicted 
level. This background level implies approximately 0.2 
background sources within r^ ■ This in turn suggests that 
essentially all of the detected sources within r^ belong to 
the cluster. We note that we have not made a correction 
for any incompleteness due to fainter sources being lost 
in the vicinity of the very bright source A. 

Given the small source sample, we have carried out 
one-parameter fits by finding the value of q which max- 
imizes the likelihood. As in our previous work, we use 
bootstrap resampling to estimate the best fit parameter 
uncertainties. The parameters r c< x and ax are expressed 
as functions of q, and So is determined by normalization. 
The resulting maximum likelihood value of the mass ra- 
tio is q = 1.17±0.28. The corresponding core radius and 
slope values are r C: x = 5.0 ± 0.8 and a c ,x = —2.5 ± 1.0. 
This fit to the cumulative radial source distribution is 
shown in Fig. 1101 For an assumed main-sequence turnoff 
mass (MSTO) of 0.80±0.05 M , the corresponding char- 
acteristic mass for the Chandra sources is 0.94±0.23 M Q . 
While this provides a hint that M x exceeds the MSTO 
mass, the difference is clearly not statistically significant. 
The large statistical uncertainty is likely due to the small 
sample size. 

We note that INovola fc Gebhardtl (|2006D have classi- 
fied NGC 6652 as core collapsed, with an optical core 
radius of 1.2". As we did in our previo us analysis of 
the C handra source distribution in M30 (jLugger et al.l 
I2007T ). we have fitted the radial distribution of sources 
with a pure power law for comparison with the opti- 
cal surface brightness profile. The resulting slope is 
ax = —1.30 ±0.22. This is significa ntly steeper than the 
slope of —0.57 ± 0.12 measured by INovola k Gebhardtl 
(|2006f ) for the optical cusp. In principle, comparison of 
these slopes can be used to estimate Mx- However, the 
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slope of the profile presented in their paper continuously 
steepens with increasing radius beyond the core, rather 
than maintaining a constant value for a significant run 
of radius as in other core collapsed clusters. The char- 
acteristic slope at 10" is about —1.3, which is the same 
as the slope for the Chandra sources. This suggests a q 
value of unity, putting Mj near the MSTO mass. Once 
again, there is considerable statistical uncertainty in this 
value. In sum, the spatial distribution of sources does 
not provide a definitive estimate of the source mass. 

4. DISCUSSION 

4.1. Source A 

Source A's unabsorbed X-ray luminosity is Lx (0.5- 
10 keV) = 7.0±° 7 7 x 10 35 ergs/s. This is compara- 
ble to our 2008 observation which found Source A to 
have a l uminosity of L x (0.5- 6.0 keV) = 4.4t° ;fj x 10 35 
ergs/s (|Coomber et al.l 120111 ; our 0.5-6.0 keV lumi- 
nosity is 5.4^q' 4 x 10 35 ergs/s). Observations of XB 
1832-330 previous to our 2008 observation dSidoli et al.l 
120081: iTarana et all l2007t iParmar et all 120011 ) have re- 
ported X-ray luminosities several times higher. We 
note that the RXTE bulge s can monitoring observations 
(|Swank &; Markwardtl 1200 if ) show a roughly sixfold 
drop in flux in early 2011. Our observation suggests 
that this decline in A's luminosity has not continued 
smoothly. 

The scheme for dividing LMXBs into ultracompact 
(white dwarf don ors) vs. long-per iod systems based 
on X-ray spectra ([Sidoli et al.l I2001D suggested Sourc e 
A (XB 1832-330) as ultracompact (jParmar et all 120011 ). 
w hich is ruled out by the 2.15 hour period measurement 
of lEngel et all (|2012f ). This scheme may still be useful, 
but in dividing short-period (< 2 hours) from long-period 
(>6 hours) systems. 

4.2. Source B 

The peak luminosity of Source B suggests that the sys- 
tem contains a neutron star; however, Source B's ex- 
treme variability is unusual for a LMXB. In May 2011, a 
6.5 hour observation of NGC 6652 taken with Gemini's 
GMOS-S CCD detector using g' and r' filters found a sim- 
ilar, h ighly variable lightcurve for Source B (|Engel et al.l 
120121) . No clear periodicity was found. 

One previously suggested explanation for Source B's 
behavior is a system with a high inclination angle, so 
that our view of the central X-ray object is obscured by 
the accretion disk. If this were the case, we would ex- 
pect Source B's spectrum to be hardened during periods 
of low count rate. We see marginal evidence for this as 
the hardness ratio ranges from 0.16±0.02 for the second- 
lowest count rate range to 0.23±0.02 for the second- 
highest count rate range. While the tendency is not 
clear, the Ng does appear to decrease with increasing 
count rate. However, the maximum Nh (2.7lJ'| x 10 21 
cm~ 2 ) is not sufficient to significantly decrease the coun- 
trate. An edge-on accretion disk, where flux variations 
are largely due to obscuration by changing Nh columns, 
can therefore be ruled out. 

The behavior of Source B may be better explained by 
an instability involving the propeller effect, when the 

12 http://asd.gsfc.nasa.gOv/Craig.Markwardt//galscan/mairi.html 



magnctosphcrc of the neutron star rotates faster than 
the innermo st region of the accretion di sk, pushing ma- 
terial away (|Illarionov fc Sunvaevl fl975l) . The location 
of the inner edge of the accretion disk will be deter- 
mined by the ratio of gas to magnetic pressure; as gas 
accumulates in the disk, the disk will come closer, en- 
abling accretion. This can produce cyclic mass trans- 
fer on a range of timescales, whic h have been explored 
in simulations by several t heorists dSpruit fc Taarr]|1993t 
iD'Angelo fc Spruitl l2010t iRomanova et all 120111 ). Such 
an instability has been suggested to explain ~ 1 Hz 
modulation in the tr ansient LMXB SAX J1808.4-3658 
(|Patruno et al.l I2009D. and the behavior of the 'Rapid 
Burster' ( Spruit fc Taamlll993l ). A similar phenomenon 
may be responsible for the rapid flaring in Source B. 

Source B has the properties of a very faint X-ray tran- 
sient ( VFXT) . These systems typically have neutron star 
or black hole primaries and have peak X-ray luminosities 
10 34 < Lx < 10 36 wi th quiescent Lx generally an order 
of m agnitude lower (|Wiinands et al.l 120061 iMuno et al.l 
2005]). Their time -averaged mass transfer rates are rather 
lower than expected for most evolutionary states. One 
explanation is that the majority of the inflowing gas 
is ejected, through e.g. a propeller effect. Weak and 
strong propeller regimes, corresponding to weak and 
strong outflows, have b e en studied by MHD sim ulations 
(|Romanova et al.l 120051: lUstvugova et al.l 120061 ). If the 
propeller effect is truly responsible for the characteristics 
of Source B, it may prove important to understanding the 
nature of other VFXTs. 

4.3. Source C 

The soft spectral component in Source C's spectrum 
is consistent with some magnetic CV spectra. Source 
C fits well to a simplified polar CV model, where a 
higher energy mekal component describes hard X-ray 
emission originating the accretion shock column front, 
while a lower energy blackbody component describes 
soft X-ray emission from the white dwarf's polar caps 
(|Ramsav et al.l [2004T ) . While Source C is several times 
more luminous than the suggested polar CV X10 in 47 
Tuc (|Heinke et al.l [20051) . its spectrum is comparable. 
Source C may therefore be the second polar CV iden- 
tified in a globular cluster. 

4.4. Source D 

Source D's spectrum was fit well to a neutron star at- 
mosphere model, while a mekal model proved to be an 
exceptionally poor fit. We found Source D to have low 
variability over our observation, while its inferred lumi- 
nosity is comparable to that seen in the 2008 observation. 
Source D's luminosity is simila r to other known quie scent 
LMXBs in globular clusters (jHeinke et al.ll2003bD . As 
well, it lies close to the theoretical cooling track for a neu- 
tron star atmosphere on the color-magnitude diagram. 
We therefore suggest that Source D is a NS qLMXB. 

4.5. Fainter Sources 

The brightest of the fainter sources, Source E's spec- 
trum and luminosity suggest that it may be a CV. How- 
ever, due to the low number of counts, the evidence is 
not conclusive. We were unable to adequately model 
Source F with any single component model, though a 
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neutron star atmosphere model proved an exception- 
ally poor fit. Sources H through K cannot be conclu- 
sively classified. However, most X-ray sources in glob- 
ular clust ers with similar lum inosities and X-ray colors 
arc CVs (jPoolev k Hutl 120061 ) . with smaller fractions of 
chromosphcrically active binaries and millisecond pulsars 
(MSPs); therefore we may suggest the majority of our 
fainter sources are CVs. Due to the hardness of Source 
L's spectrum and its distance from the core's half-mass 
radius, Source L is likely a background source. This is 
roughly cons istent with the 2.1 ba ckground sources ex- 
pected from iGiacconi et al.1 (|2001l ) extragalactic source 
counts (see Section 1X7)) . 

In 2010, NGC 6652 was detected by the Fermi Large 
Area Telescope, suggesting t hat the cluster hos ts a sub- 
stantial population of MSPs. lAbdo et al.l ([2010D have es- 
timated the number of MSPs in NGC 6652 to be 54t?,£; 
however, previous studies using r adio observations hav e 
placed this number much lower (iKulkami et al.lll990D . 
A MSP coincident with NGC 6652 has recently been re- 
ported; however, it remain s unclear whether it is truly 
associated with the cluster (|DeCesar et al.ll201lD . MSPs 
could be one or more of the X-ray sources in NCC 6652. 
Prev ious observations of MSPs in other globular clusters 
(e.g. IBogdanov et~aTl (|2005l[2010t) ) have shown that the 
brightest MSPs have hard spectra with luminosities typ- 
ically <10 32 erg s -1 . Thus, some of the fainter sources 
in the cluster may be MSPs. 



4.6. X-ray Source Population 

NGC 6652 is unusual in hosting two relatively bright 
X-ray sources. Calculations of the encounter rate, T, for 
NGC 6652 vary significantly, ranging from 0.2% to 2.3% 
of the Gal actic va l ue, de pe nding on whether values ar e 
taken from Harris! (| 19961) or lNovola k Gebhardtl (|2006l ). 
A full analysis of the effects of uncertainties in input pa- 
rameters on the encounter rate is underway (Bahramian 
et al., in preparation). A preliminary estimate of the 
uncertainty in NGC 6652 's encounter rate, using quoted 
measurements an d uncertainty in surfa c e brig htness and 
core radius from iNovola k Gebhardtl (120061) . and ex - 
tinction and distance measurements from IHarrisI (f2010h . 
yields an encounter rate T =1.5-4% of the total Galactic 
globular cluster interaction rate. Such a high encounter 
rate would explain NGC 6652 's containing two XRBs 
above 10 33 ergs/s. 

We present cluster parameters and X-ray source num- 
bers in three luminosity bins for several clusters in Ta- 
ble [4) ordered by central luminosity density. The ra- 
tio of the number of brighter to fainter X-ray sources 
(e.g. >10 33 to 10 32 " 33 erg s" 1 , or 10 32 ~ 33 to 10 31 " 32 erg 
s _1 ) in NGC 6652, and other clusters with very dense 
cores, is relatively high, indicating a relatively flat lu- 
mi nosity function. The se data support the argument 
by iPoolev et al.1 (|2002b| ) that clusters with higher cen- 
tral densities have flatter l u minosi ty functions, and the 
argument by iHeinke et al.l (|2003b[ ) that the differences 
above 10 31 ergs/s are driven primarily by changes within 
the CV population Where detailed optical followup 
has been done (e.g. iPooley et al.ll2002at Edmonds et al.l 
I2003al iLugger et al.ll2007t ICohn et aTTl2010f) . a majority 
of the X-ray sources that we consider (Lx (0.5-2.5 keV)> 



10 31 ergs/s) are associated with CVs and qLMXBsP 3 ! 
CVs make up the majority of the two lower bins, while 
the top bin is comprised of accreting NSs (qLMXBs and 
LMXBs). So the CVs tend to show higher luminosities 
on average in the densest clusters; the lack of numerous 
qLMXBs in NGC 6652 and M15 suggests that NSs are 
also affected. Overall, the data suggest that clusters that 
reach very high central density efficiently convert their 
numerous lower luminosity X-ray binaries into a smaller 
group of higher luminosity CVs and LMXBs. 

Core collapse may involve r ecurrent episodes of ex- 
tremely high central density (jFregeau et al.l I2003D . in 
which new binaries are formed and older binaries are 
destroyed or ejected through interactions. The highest- 
density clusters are also those clusters where mass- 
transferring binaries are most rap idly destroyed (or 
ejected) through further interactions (jVerbuntll2003| ). Bi- 
naries formed through interactions will tend to consist of 
two relatively massive stars. These will produce rela- 
tively high mass transfer rates, and thus relatively high 
Lx systems. Among CVs, magnetic cha nneling of accre- 
tion is required to enable high Lx (e.g. lEdmonds et al.l 
l2003bf ) : for qLMXBs, the quiescent Lx may reflect the 
time-averaged mass- transfer rate (jBrown et al-H^Sh F 3 ! 
As systems age, the companion is whittled away, leading 
to lower mass-transfer rates and thus fainter systems. 
In the densest systems like NGC 6652, X-ray binaries 
may be destroyed or ejected before they can dim signif- 
icantly. The extended radial distribution of X-ray bina- 
ries in NGC 6652, particularly the large offsets of sources 
A and B, indeed suggests recent ejection. 

4.7. Future Directions 

Source B's rapid X-ray and optical (jEngel et al.ll2012T ) 
flaring would be worth studying simultaneously (e.g. 
with ULTRACAM and Chandra), to clearly distin- 
guish the details of accretion. Further HST anal- 
ysis is needed to confirm (e.g. through star counts) 
the sharp central cusp and core-collapsed state of NGC 
6652 identified through s urface brightness studies by 
INovola k Gebhardtl (|2006| ). as well as to identify opti- 
cal counterparts of the fainter X-ray sources. If source A 
drops into quiescence during the Chandra era, a deeper 
study of NGC 6652 's X-ray sources could confirm the un- 
usual radial distribution and luminosity function trends 
we have suggested here. Detailed studies of the X-ray 
binary populations of other clusters of high central den- 
sity, together with continuing Monte Carlo and N-body 
simulations, may help us to understand how core collapse 
proceeds. 

WSS, COH k AB acknowledge support by NSERC and 
by an Alberta Ingenuity New Faculty Award. PWdetcct 
has been developed by scientists at Osservatorio Astro- 
nomico di Palermo G. S. Vaiana thanks to Italian CNAA 
and MURST (COFIN) grants. 



13 This agrees with IPoolev fc Hut) 12006T) who used the 0.5-6 
keV band. 

14 Though neutr ino losses from e.g. direc t URCA can reduce 
their quiescent Lx, Yakovlev fc Pet hick 2004. 
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TABLE 1 
X-ray Sources in NGC 6652 



Source 


R.A. 


Decl. 


Err 


Counts 
(0.3-7 keV) 


Luminosity 
(0.5-10 keV) 


A 


18:35:43.671 


-32:59:26.31 


0.29 


1267+. 3 g a 


7.0±° 7 7 x 10 35 


B 


18:35:44.564 


-32:59:38.49 


0.28 


4634l 6 6 |' 


- 


C 


18:35:45.752 


-32:59:23.24 


0.30 


596+24 




D 


18:35:45.651 


-32:59:26.17 


0.30 


631+J1 




E 


18:35:46.240 


-32:59:29.36 


0.33 






F 


18:35:47.743 


-32:59:19.76 


0.36 


35ti 




11 


18:35:45.473 


-32:59:25.48 


0.40 


16±S 


i.5t\i x 1031 


1 


18:35:45.818 


-32:59:35.82 


0.40 


isi! 


4.2+j-* x 10 31 


.1 


18:35:45.178 


-32:59:34.54 


0.43 


lL -3 


3.1+Vi x 10 31 


K 


18:35:44.078 


-32:59:00.54 


0.45 


9+ 4 

9 -3 


2.5+Q 8 x 10 31 


L 


18:35:43.891 


-33:00:45.63 


0.72 


4+3 


i.ittt x io 31 



Note. — Positions, relative position errors and counts obtained for each observed source in NGC 6652. Values 
for sources A-F were determined by running the CIAO tool wavdetect over 0.3-7 keV. Values for sources H-L were 
determined by runnin g the PWDetect sc ript over the same energy range. Positional errors are 95% confidence errors 
from Equation (5) of Hong ct al. (2005|). Values for errors in the counts arc indicative of la and are derived from 
ICehrels! i 19861) . Luminosities for sources B-F can be found in Table [3] Luminosities for sources H-L are based upon 
the number of counts (see text) and are listed in erg s . a From readout streak. 



TABLE 2 
Count rates for Source B 





Counts 




N H 


kT 


Lx 




Count Rate 


0.3-1 keV 0.3-7 keV 


Hardness 


(10 20 cm" 2 ) 


keV 


(0.5-10 keV) 


xS/dof 



< 0.05 


71 


424 


0.17±0.02 


23 + " 


>33 


5.8+Q7 x 10 33 


1.05/10 


0.05 - 0.075 


67 


410 


0.16±0.02 




>5 


1.0+° :2 x 10 34 


0.89/9 


0.075 - 0.1 


109 


679 


0.16±0.02 


2ll 9 s 


16+ 36 


lA+° i x 10 34 


0.80/18 


0.1 - 0.15 


135 


694 


0.20±0.02 


13+? 


i2ir 


x 10 34 


0.53/18 


0.15 - 0.2 


147 


720 


0.20±0.02 


wis 


7+8 
'-3 


2.4+0' 3 x 10 34 


0.48/19 


0.2 - 0.27 


107 


507 


0.21±0.02 


io±| 


>7 


ZX+H x 10 34 


0.84/12 


0.27 - 0.35 


130 


570 


0.23±0.02 


<9 


>9 


4.7l° : l x 10 34 


1.35/14 


> 0.35 


106 


491 


0.22±0.02 


<13 


131? 1 


6.5+iJ x 10 34 


0.36/12 



Note. — Number of counts observed over indicated energy ranges corresponding to particular count rate ranges. Count rate ranges 
are listed in counts/s. The hardness ratio is quantified by examining the fraction of counts in the low en ergy range (0.3 -1 keV) against 
a larger energy range (0.3-7 keV). Errors in the hardness ratio are indicative of la and arc derived from Gchrcls (1986). The hydrogen 
column density (Njj), plasma temperature (kT) and luminosity are taken from fitting the spectra over each range to an absorbed hot 
X-ray plasma emission model (XSPEC model phabs * mekal). Luminosity is in erg s . The five ranges with the highest count rate 
include a pileup component. 
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TABLE 3 
Spectral Fits 



Source 


Model 


N H 
(10 20 cm" 2 ) 


r 


kT 
(keV) 


Lx 
(0.5-10 keV) 


X^/dof 


B 
B 


POW 
MEKAL 


nil 
iot! 


1 0+0.1 

L - a -o.i 


o 4 +28 


i.7toi x i° 34 

i.6to;i x 10 34 


1.20/39 
1.13/39 


C 

c 
c 


POW 
MEKAL+MEKAL 

and A/TTPT*" A T 
£ IVI I2j x\_ i-V _L/ 

MEKAL+BBODY 
BBODY 


(5) 
<32 

<19 


1.6 


>17 

n ,o+fl. 03 
°- 18 -0.03 
>18 

07+ 01 

u ' u '-0.02 


1.3xl0 33 
1.4t°' 2 x 10 33 

9 +0.'8 v 1n 32 

1.4^0 2 x 10 33 
2.3±f;g x 10 32 


4.18/14 
1.06/11 

1.05/11 


D 
D 
D 


POW 

"ft jTTIT^ A T 

MEKAL 
NSATMOS 


42+ 14 

(5) 
q+2 
M -2 


5 1+O.8 
0.1_ 1 2 


1.3 

n iic- +0.002 


4.1^ x 10 33 

8.3x 10 32 
9.0+ 4 ;? x 10 32 


0.77/14 
3.52/15 
0.65/15 


E 
E 


POW 
MEKAL 


cq+48 
do -36 
04+35 
°*— 26 


2 2+ 1 ' 1 
z,z -0.8 


>1.8 


3.8 +3 i x 10 32 
2.9t° i x 1° 32 


0.62/6 
0.78/6 


F 
F 
F 


POW 
MEKAL 
NSATMOS 


<35 
<82 
(5) 


1 7+ 1 - 1 


>1.9 

8.5xl0- 2 


1.1+05 X 1° 32 

lA±°i x 10 32 
1.3xl0 32 


1.03/1 
1.11/1 
6.54/2 



Note. — Spectral fits for X-ray sources in NGC 6652. Sources B through E arc fit from 0.5-10 keV; however, due to a low number 
of counts (42) Source F is fit from 0-10 keV. All models included absorption (XSPEC model phabs). For poorly fit models (x 2 > 
2) Nh was kept at the cluster value and errors were not calculated. The luminosity stated is for the unabsorbed X-ray luminosity 
in erg s . Two component models are continued on a second line and give the X-ray luminosity separated by components. The 
total unabsorbed 0.5-10 keV luminosities for Source C from the MEKAL+MEKAL model and the MEKAL+BBODY model 
are respectively 1.8tg'2 x 10 33 er S s_ 1 an d l-^to 2 x 10 33 cr S s — 1 . The final column lists the reduced chi-squared and the degrees 
of freedom for each model. 
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TABLE 4 

Parameters and X-ray Sources of Selected Clusters 



Cluster 


Distance 


Pa 


r c 


r 


Number of X-ray Sources 


References 




(kpc) 


(L Q pc-3) 


(pc) 


(% Gal.) 


(> 10 33 ) (lO 32 " 33 ) (i0 31 ~ 32 ) 





47 Tuc 


4.5 


7.8 x 10 4 


0.46 


2.9 


1 


3 


27 


1 


M28 


5.5 


1.3 x 10 5 


0.26 


2.0 


1 


2 


14 


2.3 


M80 


10.0 


1.4 x 10 5 


0.22 


1.6 





1 


12 


1 


NGC 6440 


8.5 


1.7 x 10 5 


0.35 


3.6 


1 


10 


>13 


5.3 


NGC 6388 


9.9 


2.3 x 10 s 


0.21 


3.1 





13 


27 


6 


NGC 6752 


4.0 


2.3 x 10 s 


0.13 


1.2 





1 


8 


7 


M30 


8.1 


3.6 x 10 5 


0.063 


0.6 


1 





3 


8 


NGC 6397 


2.3 


4.1 x 10 5 


0.041 


0.3 





2 


7 


9 


NGC 6652 


10.0 


1.0 x 10 6 


0.058 


2.3 


2 


3 


■1 


10 


M15 


10.4 


3.5 x 10 6 


0.050 


6.7 


3 


2 


>2 


11,12 



Note. — Parameters and X-ray source numbers from clusters with numerous published X-ray sources above Lx = 10 31 ergs/s 
(Terzan 5 omitted due to uncertainty in parameters). The numbers of X-ray sources include only sources within the cluster half-mass 
radius, using luminosities in the 0.5—2.5 keV range from the deepest observation (or average published luminosities for M15 sources). M15 
and NGC 6440 are incomplete in the lowest bi n. Central luminosity dens ities (for the core-collapsed, we average wi thin the core radius) 
and core radii are calculated using values from Noyola & Gcbhardt (2006) where given, while other values arc from Harris (2010). Close 
encounter rates T oc , 5 r 2 are given as a percentage of the total Galactic globular cluster system rate. References.- (1) IHeinke et al.l 
(2005); (2) Becker ct al. (2003); (3) Hciiike ct al. (20 03bD; (A) IHeinke et al.l (1200331 1: (5) IPoolev eTall l|2002Hl: (6) [M axwell c t~aTl pH2T ); 
(7) [Po oley et al. (2002a); (8) Lueecr ct al. (2007): (9) Grindlav et al. (20013); (10) this work; (11) [Hannikaincn ct al. (2005); (12) 
IHeinke et al.l ( 120091)" 
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